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As a quantitative measure of localization, the inverse participation ratio of low lying Dirac eigenmodes and 
topological charge density is calculated on quenched lattices over a wide range of lattice spacings and volumes. 
Since different topological objects (instantons, vortices, monopoles, and artifacts) have different co-dimension, 
scaling analysis provides information on the amount of each present and their correlation with the localization of 
C$ ' low lying eigenmodes. 



1. INTRODUCTION 

With modern computational power has come 
the ability to examine the low lying eigenvectors 
of the Dirac operator and hence their spatial cor- 
relation with instantons and other related objects 
thought to be involved in chiral symmetry break- 
ing and confinement |H2j . While these studies 
focused primarily on the local relationship be- 
tween instantons and low-lying Dirac eigenmodes 
(LDEs), other models of confinement and chiral 
symmetry breaking involving objects of lower co- 
dimension are popular, based on monopoles, vor- 
tices, and hybrid objects 0. Presumably these 
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objects would have a rather different effect on 
the LDEs than 4-dimensional instantons, due to 
their different co-dimension. Furthermore, a re- 
cent study 0] has suggested a dense layered 3- 
dimcnsional structure to the LDEs. 

One difficulty is the quantitative characteriza- 
tion of localization of the LDEs or related quan- 
tities such as the topological charge density. In 
localization of the LDEs was studied using 
the inverse participation ratio (IPR) which yields 
a number characterizing the localization of an 
eigenmode. 

By studying the scaling dimension of the IPR, 
we can find the co-dimension of the structures 
which localize the LDEs, thus giving some insight 
as to the possible confining objects and mecha- 
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2. INVERSE PARTICIPATION RATIO 
(IPR) 

The IPR of a normalized field pt{x) is defined 



i = nJ2p1(x) 



(1) 



where N is the number of lattice sites x. Here 
we use pi(x) — iplipi(x) and ipi(x) is the i-th, 
normalized (^2 x Pi{x) — 1), lowest eigenvector of 
the Dirac operator. 

With this definition, / characterizes the inverse 
"fraction" of sites contributing significantly to the 
support of p(x) (we now drop the subscript i). A 
simple calculation shows that the IPR takes the 
following values for these simple situations: 

Unlocalizcd : p(x) = const. 1=1 
S — function : p(x) = S(x ) I = N 
localized on fraction / of sites : I = 1/f 

Suppose that the objects responsible for con- 
finement, or indeed any physics governing the 
lowest Dirac eigenmodes, localize the LDEs. As 
the lattice spacing is reduced, the fraction of sites 
contributing to the IPR scales as a d /a 4 . Thus the 
IPR indicates the co-dimension of these objects: 
d = 4 for instantons, d = 3 for monopoles, and 
d = 2 for vortices. Gauge dislocations should 
contribute as d = objects, however their den- 
sity diverges as a~ 4 so that they should give a 
^constant contribution: a°/(a 4 a~ 4 )). 

Since the IPR ~ 1//, if we reduce the lattice 
spacing at fixed physical volume, we have 

4-d 



at fixed volume : / ■ 



(2) 



On the other hand, increasing the volume at fixed 
lattice spacing includes proportionately more of 
the confining objects, whatever their dimension. 
Thus we expect the IPR to remain constant, 



L — > oo at fixed a : I ~ constant 



3. RESULTS 



(3) 



We have explored these two regimes using 
quenched lattices generated with the tadpole im- 



proved Symanzik gauge action, and the parame- 
ter set shown in Table 1. On each lattice we com- 
puted the lowest eight eigenvectors of the Asqtad 
Dirac matrix. 

Table 1 

Lattices analyzed 



a 


L 


vol 


P 


no. configs. 


a^O: 


0.20 fm 
0.15 
0.12 
0.095 


12 
16 
20 
24 


(2.4 fm) 4 
(2.3 fm) 4 


7.56 
7.847 
8.109 
8.456 


100 
97 
93 
118 


L — > oo: 


0.12 fm 


12 
16 
20 
24 


(1.4 fm) 4 
(1.9 fm) 4 
(2.4 fm) 4 
(2.9 fm) 4 


8.109 


100 
100 
93 
100 



We see clear evidence for lower dimensional 
scaling as a — > with co-dimension between 2 
and 3. Figure 1 shows both the distribution of 
the IPRs and the scaling of the averages. The 
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Figure 1. Scaling of the average IPR as a — > 

data points are fit to c\/a n + C2, where c\ and 
C2 are constants and n = 1,2, 3. The reduced chi 
squared values 2 for the fits are 6.7, 3.2, and 45 
for n =1, 2, and 3, respectively. 

2 The error bars shown here are corrected from those (much 
larger) shown at the conference. 
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In figure 2 we show the behaviour of the IPR 
as we increase the volume at fixed lattice spacing; 
it is rather unaffected, as we expected above. 

IPR scaling as "L -> infinity" 
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a=0.15 L=16 
a=0.12 L=20 
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HYP Sweeps 



Figure 3. Avereage topological IPR vs. number 
of HYP smearing sweeps, as a — » 
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Figure 2. Scaling of the average IPR as L — > oo 
with a — 0.12 fm 

4. TOPOLOGICAL CHARGE 

We can also investigate the localization of the 
topological charge density by computing the IPR 
from q(x) = F^F^, where we have normalized 
J2 X \l( x )\ = 1- We have computed q(x) by suc- 
cessive HYP smearing sweeps on the a — > 
series, and show the results in Figure 3 (note that 
only 5 HYP smearing steps were performed on the 
a = 0.12 fm L = 20 lattice set). While this plot 
does not show us new information on the local- 
ization of the LDEs, it is nonetheless instructive. 

First, we see that all lattices without smoothing 
have an IPR = tt/2. This is the value expected if 
the field is a gaussian fluctuation at each site, re- 
gardless of its width. We further see the approach 
to a stable localization of topological charge ver- 
sus HYP smearing as the lattice spacing is de- 
creased (at fixed volume). We note that <IPR> 
is not large, meaning that q{x) is not strongly lo- 
calized. Also, it increases as a — > as for the 
LDEs. 

5. CONCLUSIONS 

The main result of this study is the indica- 
tion of a localization of the low-lying Dirac eigen- 
modes on surfaces of co-dimension between 2 and 



3, qualitatively supporting the center vortex or 
monopole pictures of confinement. Note, how- 
ever, that the singularities of thin objects (vor- 
tices or monopoles) are expected to be smoothed 
out by the QCD interactions and become thick, 
with a size ~ 1/Aqcd- Thick objects fill a fixed 
fraction of space, not a divergent one. The indi- 
cation we have, via the divergence of <IPR> as 
a — > 0, of localization on singular manifolds, is 
remarkable, whatever these manifolds are. 
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